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Freflico 


This  thesis  is  the  result  of  an  invest i;'at Ion  of  auad- 
ratur-e  clock  modulation  (iJCM)  versus  binary  and  quadr'ature 
phase  shift  key  in, q (BFSK  and  CJFSK)  in  the  presence  of  inter- 
symbol  interference  (ISI).  QCM  is  a modification  of  normal 
BPSK  in  which  the  alternate  data  bits  are  transmitted  on  an 
ortho/’:onal  carrier.  QCM,  when  compared  to  BPSK,  was  found 
to  improve  the  data  performance  of  a communication  system 
which  is  limited  by  intersymbol  interference.  QCM,  when 
compared  to  QI’SK,  was  found  to  lie  sli,‘'htly  inferior  in  data 
performance.  The  carrier  trackin;"  characteristics  of  QCM 
and  BPSK  wer*’’  also  i r.ves t i,Tated  in  the  presence  of  131,  and 
QCM  was  found  to  have  the  same  carrier  track  in, q character- 
istics as  PPCK.  In  addition,  the  ener/’-y  penalty  for  track- 
ing: a QFSK  sirnal  versus  a binary  modulated  sii'.nal  was  ob- 
tained under  no  I3I  cot\ditions.  This  peiialty,  tc, ''ether  with 
the  almost  equal  data  performance  of  QCM  and  QPSK,  suz-^i'^est 
that  QCM  could  possibly  be  used  as  a back-up  systt'm  to  QPSK. 

I would  like  to  thank  my  committee  members.  Captains 
Grerr.  Vauphn  and  T.R.  Hadley,  for  their  valuable  .assistance 
and  advise  rendered  to  me  durin't  my  work  on  this  thesis.  1 
would  particularly  like  to  thank  my  thesis  avlvisor.  Captain 
Stanley  Robinson,  for  the  enthusiasm,  cu id.r-t\ce , anil  valuable 
Insi^’iht  th.at  he  provided. 

None  of  this  thesis  could  have  been  accomplished  if  t\et 
for  the  love  and  consideration  that  I reeeived  fi  om  my  fam- 
ily. First,  to  my  sons  Allen  and  Pavid,  who  kept  the  TV 
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shift  keyinr  (PPSK)  technique  in  which  alternate  bits  are 
transmitted  orthoronally . Compared  to  normal  PPSK  modula- 
tion, jCM  improves  the  performance  of  a communication  sys- 
tem which  is  limited  by  intersymbol  interference  (ISI). 
Compared  to  quadrature  phase  shift  keyinr  (CJFCK)  for  the 
same  information  rates,  QCM  is  sli^tbtly  inferior  in  data 
detection  performance.  The  average  carrier  tracking  cap- 
ability of  a CCr>'.  sc.hi-.mo  is  equal  to  that  of  PFSK  for  the 
same  closed-loop  bandwidth  and  time  bandwidth  product  (BT) 
of  the  channel  filter  at  a specified  signal  energy  to  noise 
ratio  (Eb/'>Jo). 

The  data  detectio.n  performance  for  QCV.,  PFSK,  and  JPSK 
is  analyzed  by  comparing  the  one-shot  probability  of  error 
conditioned  on  a phase  error  as  a function  of  Eb/No  and  BT 
for  a specifiea  cha.nnel  filter.  Carrier  track  itig  perform- 
ance for  3CK  and  BFSK  is  analyzed  by  obtaining  an  average 
phase  error  variance  for  the  linear  model  of  a Costas  loop. 
Tn  computing  both  the  probability  of  error  and  phase  error 
variance,  the  intersymbol  interference  is  modeled  from  a 
truncated  data  sequence.  In  addition,  the  bit  energy  pen- 
alty for  tracking  a QFSK  signal  versus  a BFSK  or  .iCV.  signal 
with  no  ISI  is  examined. 

CCM  offers  an  advantage  in  data  detection  performance 
over  BFSK  for  low  B?  products,  and  an  apparent  advantage 
over  QFSK  in  better  carrier  tracking  in  the  presence  of  a 
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QrAD.-UTl'R’r:  CLCJK  r.Or'.’l.A TIC:;  VHKJ'.O 
BINARY  AND  r AD;UrVHD  DRADr;  DHir:  K-YIND 
IN  THb:  DRDSKNCD  OR  INTDRoY: D'OL  I Nl’RRr YKRNCD 

I.  IntroA’.ic'tii'':'. 


The  probaMlity  of  orroi'  of  u Aata  corinainioat  ior:  ayaton 
can  bo  made  arbitrarily  a.rall  .'Ivon  ot'.ourb.  airriai  or.or-y  or 
sif^'-al  duration.  '.'owover , ir'.  the  real  world,  tb.e  oo.-.r  tra  into 
of  limited  aienai  enerriea,  ittoroarin.y  i:'.  format  i or.  rater,  Luid 
finite  bandwidtlir,  pl.iee  irerer  limit:;  on  the  reT'forma.nee  tt.ai 
can  be  obtaiti.ed.  In  aatoliite  cor:i.mur. ica t lone  , data  rate;; 
must  be  hieh  because  of  the  larme  volu.me  of  i-.'.form.it  ion  ";’'..it 
is  trailer,  i t ted  . Also,  the  fre-]ueney  a 1 loe;\.  t ions  for  a svi'- 
tem  place  a definite  limit  on  the  ch.uir.el  bar.dwidtb.s  t.hai 
can  be  used.  This  two-pronrod  constraint  of  increas  ine  d.ita 
rates  and  fixed  channel  b.indwidths  becomes  .irtn i f e s t in  .i 
phenomen.a  kiiown  as  ir. tersyr.bol  interference  There- 

fore,  alternative  s i-rial  1 in.:  schemes  .mu;;t  be  sou.n  t that 
yield  the  best  system  performance  for  the  lovsest  cost. 

Intersymbol  interference  is  caused  by  imissin,'  a dat.i  s-'- 
quence  through  a filter  of  finite  bandwidth.  This  effect 
can  usually  be  r.eeleoted  wb.en  tb.e  bandwidth  of  the  filter 
is  much  laryer  than  the  dat.\  rate.  Howe\er,  as  the  data 
rate  appro.aches  the  bandwidth  of  the  filter,  the  effects  of 
ISI  become  more  pronounced.  The  de, -'rad  ins  effects  from  IS! 
are  twofold.  First,  the  power  is  reduced  in  a y.iver.  datum 
epoch  by  the  smearin.t  and  distortion  caused  by  the  filter. 
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If  not  compensated  with  an  increase  in  transmitter  power, 
the  resulting  performance  of  the  system  will  suffer.  Sec- 
ond, reinforcement  or  cancellation  of  the  desired  decision 
variable  from  the  other  data  pulses  occurs  on  a random  basis. 
The  net  result  of  this  random  process  is  reduced  system  per- 
formance . 


Various  methods  have  been  employed  in  an  attempt  to  re- 
duce the  effect  of  intersymbol  interference.  Linear  equal- 
izers can  reduce  intersymbol  interference;  however,  some 
additional  noise  is  introduced  that  reduces  the  total  ef- 
fective signal- to -noise  ratio  of  the  channel  below  that  of 
an  ideally  wide  channel  (Ref  1).  Nonlinear  decision  feed- 
back equalizers  introduce  less  noise;  however,  they  only 
consider  past  data  bits.  Other  techniques,  such  as  partial 
response  and  coding  schemes,  have  also  been  used  to  reduce 
the  effects  of  I3I. 

Another  iriethod,  quadrature  clock  modulation  (QCM),  in 
which  the  alternate  data  bits  of  a binary  phase  shift  keyed 
(BPSK)  system  are  transmitted  orthogonally,  has  been  sug- 
gested as  a means  to  reduce  the  effects  of  ISI  (Ref  2). 

A comparison  of  QCM  to  BPSK  has  been  made  (Ref  3)* 
However,  because  of  the  increasing  emphasis  on  multiphase 
systems  and  since  QCM  requires  the  additional  complexity  of 
a quadrature  channel,  QCM  should  be  compared  to  quadrature 
phase  shift  keying  (QPSK)  (Ref  ^033)  ♦ 


Background 


The  Defense  Communication  Agency  (DCA),  the  directorate 
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of  tho  l'(' fi'tnu'  o;\ it' I I i t.f'  v'ommun\o;itvon  (l):iCo),  h;i 


propo:;«'il  ^I'oK  anil  I'l'oK  aa  tho  modulation  loohnlmioa  for  l);'>Co 
Fhar.o  111.  Th«'  r.yrtom  will  \itlll;.o  anrfioo  aoonr.tlo  wavi' 
(SAW)  fllti'ra  with  a oomhlnovl  rooolvi'r  and  transmit  lor  I'nnd- 
width  tlmo  prodno  t of  1.0  lK*'f  *')•  Tin*  Impulr.o  oharao  t«M' I r- 
tica  of  tho  SAW  flltorr.  Ir.  rouphly  I'lpilvalont  to  a flvo-polo 
Puttorworth  flltor.  Nlolaot\  modolod  thor.o  r.yr. ti'm  flllo!".?  ar. 
a alx-polo,  0.1  dh  rlpplo  Olvd'yahov  I'andpar.s  flltor  aiul  o\al 
natod  tho  por  formatiO>'  of  'd'SK  vorr.ur  vjOM  for  !'T  produo  tr.  of 
0.*',  1,  aiAil  l.s.  Ih'oa'uu'  Id'SK  la  oaaoi\tlally  tho  haok-up 
modulatlot\  r.ohomo  for  USCS-lll,  ah.oulil  alr.o  ho  oompai'ovl 

to  Ql'SK  to  aoo  If  at»y  lmprovomot\ t oat\  ho  mado  In  roduoinr, 
tho  of  foot;’,  of  ISl. 


lt\  nddltloi\,  tfio  oarrlo!-  traolvli^r,  por formunoo  for  a viOM 
syatom  haa  i\ot  l'<'Oi\  i tn-oa  1 1 ya  tml . It  la  i\ot  kt'iowr\  If  thla 
will  Introduoo  additloiial  prohlomr.  boyond  tho  oarrlor  traok- 
In^’;  dlfflcultloa  oi\ooui\torod  In  a !U'SK  or  vjrSK  ayatom. 


Frobl om  Slat  omon  t 

Tho  pui'poao  of  thla.  thoa.la  la  to  d»'toT*mlno  tho  I'olatlvi' 
porformanoi'  of  id'SK , OdM,  and  .Jl'SK  ut\dor  t lmo  and  bat\d\sldth 
cona  tra  In  t r.  th.vt  produoo  a I y.ti  I f io.vn  t amounta  of  In  ti'rayml  o 1 
intor  foronoo  i and  to  I tivoa  t I pa  to  tho  oarrloi-  traoklnp  por 
formanoo  of  lidM  vo;-a»ia  bl’SK  and  lii'SK . 

ProMom  lji’j;_\_a 

Tho  proolom  analya.la  la  dlvldod  Into  tht'oo  poi\oral  aroar- 
Flrr.t,  thi'  r.yr.tom  doaorlpf  Ion  of  tl\o  a.lpn.il  aot  and  data 
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detection  and  carrier  track inr.  component;5  I’or  the  thi-ee  mod- 
ulation schemes  are  discussed.  Second,  the  various  view- 
points of  analyzinr  intersymbol  inter t'erence  are  examined. 
Third,  the  scope  and  limitations  of  the  problem  resolution 
are  stated. 

System  Descr  irt  ion  > The  si/m.al  sets  for  the  three  mod- 
ulation schemes  are  shown  in  Fipure  1.  The  si,ynal  set  for 
QPSK  call  be  thought  of  as  two  independent  binary  channels. 

For  this  reason,  t;rSK  has  half  the  data  or  symbol  bandwidth 
of  BPSK  for  the  same  information  r.ate. 

Quadrature  clock  modulation  is  a binary  modulated  trans- 
mission scheme  in  which  the  alternate  data  symbols  are  ti\uis- 
mitted  on  an  ortho,;onal  signal.  Thus,  QCM  has  the  same  data 
bandwidth  as  DP3K  for  the  same  information  rate.  In  tliis 
thesis,  a BT  of  2a  for  QPSK  corresponds  to  a BT  of  a for 
BPSK  and  QCM  at  the  same  information  rates  and  usin,'t  the 
same  channel  filters. 

A modulation  scheme  for  QCM  is  shown  in  Figure  2.  The 
logical  data  sequence,  the  I sipnal,  is  transmitted  alter- 
nately on  the  sine  or  cosine  channel.  For  this  thesis,  the 
sine  channel  corresiionds  to  the  odd  data  pulses  and  the  co- 
sine channel  corresponds  to  the  even  data  pulses.  Because 
the  data  is  transmitted  orthcr.onally , the  effects  of  inter- 
symbol interference  are  reduced.  In  particular,  the  ad- 
jacent interference  of  the  previous  and  past  data  pulses 
are  transmitted  on  an  orthogonal  carrier  and  are  not  ap- 
parent to  the  data  detector  if  perfect  phase  coherence  can 
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be  maintained.  Even  if  perfect  phaac  coherence  cannot  be 
maintained,  the  effectr.  of  I3I  can  r.tili  be  a ipTii  f ican  tiy 
reduced  for  phase  errors  of  less  than  n/h  radians. 

The  data  detector  for  those  modulation  scnomes  is  shown 
in  FifTiirc  3*  Foi'  BPSK  and  QCM  the  output  of  the  VCO  is  the 
reference  si^mal.  In  the  case  of  QCM,  the  reference  si^-jial 
is  shifted  by  tt/2  radians  with  every  datum  pulse.  This 
clockin.y  of  the  reference  sirnal  causes  the  data  to  always 
appear  on  the  cosine  le,'’  of  tiie  Costas  loop. 

Since  BFSK  and  QCM  arc  binary  modulated,  the  carrier  can 
be  tracked  as  shown  in  FI, ''cure  3*  QFSK,  however,  with 

its  quadrature  modulation  set,  a reference  sii'Sial  external 
from  the  data  detector  must  be  provided  to  resolve  tiie  mod- 
ulation ambirnities  and  recover  the  carrier.  A Costas  loop 
that  tracks  a quadrature  sifTial  (N  - ^ Costas  loop)  is  sho.vn 
in  Fiivure  4. 

Intersvmbol  Trrter  fere  nee . The  effi'cts  of  intersymbol 
interference  otr  data  performance  and  cai'rier  track  in, r ia  the 
main  subject  of  this  thesis.  The  problem  of  system  evalua- 
tion is  how  to  calculate  the  probability  of  error  or  carrier 
tracking;  parameters  so  as  to  include  the  effects  of  inter- 
symbol interference. 

For  the  data  detection  problem,  the  major  efforts  have 
been  in  two  areas.  First,  various  authors  have  attempted 
to  find  density  functions  or  bounds  on  the  probability  of 
error  which  include  the  effects  of  intersvmbol  interference 
(Refs  6 , 7 , 8 , , 10 ) . Four  types  of  ’’ovinds  have  been  suf^qes  tod 
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(Ref  9).  The  worst  case  bound  assumes  that  all  intorferers 
add  in  phase?  therefore,  this  bound  is  too  pessimistic  to 
represent  a data  detector.  The  second  is  the  Chernoff  bound 
(Refs  ?.?)•  These  bounds  are  difficult  to  obtain  when  the 
number  of  s i, ‘'nail in?’;  phases  is  rreater  than  four.  Also,  the 
Chernoff  bounds  tend  to  be  a fairly  loose  approximation  and 
become  more  inaccurate  as  the  number  of  poles  in  the  filter 
increases  (Ref  Another  metliod  is  by  a pair  of  con- 

verf'in?'  bounds  where  a lower  and  an  upper  bound  are  found 
in  terms  of  mar/rinal  distributions  (Ref  ^) . These  paired 
bounds  are  difficult  to  obtain  as  the  number  of  s i?:nall in?' 
phases  increases.  Also,  they  do  not  model  or  represent  the 
errors  due  to  crosstalk  between  the  in-phase  and  quadrature 
channels  that  appear  in  QPSK.  The  fourth  type  of  bo’und  as- 
sumes that  intersymbol  interference  can  be  treated  as  addi- 
tive Gaussian  noise  with  the  same  power.  However,  work  by 
Shimbo,  Fanp,  and  Celebiler  discount  this  last  method.  Den- 
sity functions,  approx imat in?'  intersymbol  interference,  have 
also  been  f;enerated  (Ref  10).  The  density  function  method 
appears  to  have  merit?  however,  based  on  observations  of 
different  pulse  responses,  this  author  questions  the  par- 
ticular series  representation  of  the  inters^mbol  interfer- 
ence. It  would  appear  that  the  coefficients  in  the  scries 
would  have  to  be  matched  to  the  specific  filter  be  inn  modeled. 

The  second  approach  that  has  been  suppested  is  to  numer- 
ically calculate  a finite  number  of  i.ntorference  terms  (Ref 
11 ).  The  major  limitinfr  factor  in  this  method  is  the  number 
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of  computa  t ioru:  iiK'ri'ai'.oa  oxpot'.ont  ially  aa  the  level  of  mod- 
ulutlot\  i'i.-nala  ita^reaae.  The  finite  calculation  nu'thod  io 
alao  an  approx  i:iia  t ion  to  an  exact  rol\itio:\  which  contains 
an  infit\ite  tuanl't'!*  of  terma.  Fo!-  hPilK  and  ^I'SK  the  co;nhit\- 
atorial  por.r.  ih  il  i ties  are  within  reason,  and  finite  calcu- 
lation will  ht'  ur.t'd  foi’  this  thesis.  Also,  this  method  of 
numerical  r.  imulatior^  letids  itself  to  easily  evaluate  an  ex- 
press iot\  for  the  phase  error  varia!\ce  of  tlie  '.'1'Sa  or 
carrier  track  ins  loci'  in.  th.e  pt'esenct'  of  ii'. tersymhol  ir.ter- 
ferer.ce  . 

With  this  method  it  is  necessai'y  to  determine  how  man.y 
bits  of  a truncated  sei]ue:\ce  to  inclvuie.  Jones  has  found, 
except  for  r.mall  PT  products,  r.a  t is  f.ic  tony  rt’sul  ts  can  'ee 
obtained  with  the  adjaceiit  intt'rsynbol  in. terference  terms 
(Ref  V'o!-  thi^•.  thesis,  the  d.it.i  stream  will  be  trut'.- 

cated  so  an.  to  include  inter’symbol  intt'i'ferers  that  are  with- 
in a specified  power  threr.b.old  relative  to  the  major  datum 


pulse . 

Problem  b'olutlor'. . I'ht'  major  objective;;  tb.at  will  be 
resolved  by  this  thesis  a’.'et  (1)  Jompai'e  I'PJK  to  dOM  for 
various  RT  products,  and  s i.nial- to-no  ise  ratios.  (f)  Com- 
pare dPSK  to  dCy.  for  the  sanu'  information  rates,  filtei' 
bandwldttu;  at'.d  transmitter  powers.  (li  Obtain  express  iotu; 
for  phast'  error  variance  of  a carrier  track  in,’;  loop  for 
PPSK  and  d'-l"!  in  the  prest'nce  of  IJI.  (a'  Compare  the  enersy 
reqviiremet\tn.  for  the  carrier  track  iny  of  JCM  and  RFSK  vennus. 
QPSK. 


11 


Fundamental  Assumptions 


Assumptions  must  be  made  to  narrow  the  problem  to  one 
which  can  be  easily  solved.  Those  assumptions  .'lobal  to 
the  entire  system  are  riven  here;  assumptions  needed  to  de- 
velop a specific  equation  or  concept  cf  the  mathematical 
models  will  be  stated  as  required  in  Chapter  II. 

Because  data  performance  and  carrier  tracking  are  in- 
terrelated, an  exact  solution  would  require  a simultaneous 
nonlinear  analysis  of  the  problem.  However,  for  practical 
applications,  the  carrier  tracking  performance  can  be  sep- 
arated from  the  data  detection  problem.  The  basic  assump- 
tion to  this  approach  is  that  the  carrier  tracking  device 
is  initially  in  lock,  at  some  offset  phase  error,  and  re- 
mains constant  over  the  range  cf  specified  s ignal- to-no ^se 
ratios.  Another  way  of  stating  this  is  that  the  carri'^r 
tracking  device  continues  to  track  for  a few  d3  of  signal 
to  noise  de^rredation  beyond  which  the  data  detector  fails. 

The  data  detector  of  Figure  3 is  assumed  linear  so  su- 
perposition is  valid  and  a baseband  model  can  be  analyzed. 
The  phase  errors  are  assumed  to  be  small  so  the  carrier 
tracking  loop  can  be  replaced  by  a linear  model  as  shewn  in 
Figure  5*  Also,  exact  timing  information  is  assumed  for 
proper  operation  of  the  integrat-j  and  dump  circuits  and  for 
the  Tr/2  phase  shift  of  the  VCO  output  for  QCM  detection. 

The  noise  process  added  to  the  signals  prior  to  detec- 
tion is  a zero-mean,  stationary,  Gaussian  random  process 
with  a two-sided  power  spectral  density  cf  No/2  (watts/'Hz). 
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The  final  assumptions  concern  the  data.  Each  bit  of 


continuous  data  is  assumed  to  be  statistically  independent 
and  identically  distributed.  The  one-shot  probability  of 
error  is  computed  assum.i.np:  an  uncoded  data  stream  for  both 
the  EPSK  and  QCK  cases.  For  QPSK,  a Gray  code  is  assumed 
so  as  to  compare  a QPSK  bit  probability  of  error  to  the 
other  cases.  Also,  any  phase  ambiguities  in  phase  tracking 
arc  assumed  to  bo  resolved  by  differential  encoding  of  the 
data. 

Research  Ke  thod 

First,  all  of  the  necessary  mathematical  models  will  be 
developed  to  answer  the  major  objectives  of  this  thesis. 
These  models  will  then  be  programmed  for  numerical  analysis. 
Case, results  for  data  detection  performances  will  be  ob- 
tained for  QCM  versus  EPSK  for  BT  products  of  0.5,  1*0,  and 
1.5  and  with  fixed  phase  errors  of  0°,  10°,  and  20°.  QCM 
and  QPSK  will  bo  compared  at  the  same  information  rate,  fil- 
ter bandwidth  and  transmitter  power  for  BT  pairs  of  (0.5, 
1.0),  {0.75 1 l-5)f  and  (1.0,  2.0).  The  carrier  tracking 
phase  error  variance  for  3PSK  and  QCM  v/ill  be  given  for  3T 
products  of  0.5i  1*0,  1.3>  and  infinity.  Finally,  the  car- 
rier tracking  penalty  for  a QPSK  system  under  no  ISI  con- 
ditions will  be  investigated. 
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filters.  The  vise  of  numerical  analysis  dictates  a d.ecisicn 
for  the  model  of  the  filter.  DCA  has  proposed  a sv'ctem  tivat 
utiliv.es  surface  acoustic  wave  filters  tlvat  have  a response 
similar  to  a five-pole  fnitterworth  filter.  However,  the  im- 
pulse response  of  these  filters,  was  not  available  for  use  in 
this  thesis.  It  Ivas  been  shown  that,  in  general,  inters.yra- 
bol  inter ferenoe  increases  with  an  Increase  in  the  tuimber  of 
poles  for  both  Chebys'nev  and  Butterwor tlv  filters  (Hef  11). 
With  the  philosophy  of  obtctinin^t  a close  upper  bound  by  in- 
creasing: the  poles  by  one,  Nielsen  chose  a six-pole  Cheby- 
shev  filter  to  model  the  cb.aivnel  filter. 

In  order  to  correlate  the  results  of  tivis  thesis  with 
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those  of  Nielsen,  the  filter  response  that  ho  chose  will 
also  be  used  here.  The  transfer  function  of  the  six-pole, 

0.1  dB  Chebyshev  filter  is 

GbO 

H(s)  ^ (1) 

s^  + b5s^  + b4s^  + b3G^  + b2s‘‘  + bis  + bO 

where  G is  the  rain  and 

bO  - 0.20713  b3  = 2.77008 

bl  = 0.90176  b4  = 2.96575 

b2  rr  2.04784  b5  = 1.71217 

and  the  frequency  is  in  radian/socond  (Ref  12t290,  Table  A. 

2). 

The  pulse  I’esponse  of  the  filter  is  obtained  bv  usinr 
two  computer  prorrans.  Factors  of  the  denominator  poly- 
nomial are  obtained  from  POLY  (a  local  polynomial  evaluation 
routine).  These  factors  are  then  used  to  evaluate  the  step- 
response  H(s)/s  with  PARTL.  PARTL  is  a local  Hcavysidc  par- 
tial fraction  expansion  and  time  response  propram.  Both  of 
these  proprams  are  briefly  discussed  in  Appendix  A.  The 
step- response  of  the  filter  with  unity  pain  is 

r(t)  = 1.0 

+ 0.94211  oxp(-0.31334t)  sin(0.773'39t  + ll6.l6S°) 

+ 0.25237  exp(-0.1l469t)  sin(1.0565t  - 19.727°) 

■f  2.32050  exp(-0.42806t)  sin(0.28310t  + 229.341°)  (2) 

The  pulse  response  is  obtained  ‘’rom  p{t)  - A(r(t)  - r(t-T)), 
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where  A is  the  normalized  voltage  narnitude  of  the  response 
and  T is  the  pulse  time  duration;  A will  be  varied  to  ob- 
tain appropriate  simal  to  noise  ratios  for  the  plots. 

Equation  2 is  modified  to  obtain  a step-response  in 
terms  of  3T  products.  Applying  the  frequency  scalinp  facxor 
in  the  form  t = DTZ-rrt',  where  t’  is  the  new  time  variable, 
yields  an  expression  in  terms  of  BT  products.  The  BT  prod- 
uct is  defined  to  be  one  when  the  one-sided  bandwidth  is 
one  Hz  and  the  pulse  duration  is  one  second. 

The  pulse-responses  for  three  BT  products  for  this  fil- 
ter are  shown  in  Fiiture  6.  This  fipure  illustrates  the  in- 
creased distortion  or  spreading  caused  with  lower  BT  pro- 
ducts. Also  apparent  is  the  celay  that  must  be  applied  to 
the  intonate  and  dump  circuits  to  maximize  their  outputs. 
This  time  or  proup  delay  can  be  obtained  in  a nvoiber  cf  ways. 
Jones  evaluated  the  derivative  of  the  phase  function  of  the 
filter  (Ref  11 1 126,  Eq  4l ) . However,  Hansell  supfiests  the 
use  of  a finite  difference  method  whenever  a baseband  model 
is  employed  (Ref  13).  Both  methods  are  suboptimum  in  ob- 
taininf  the  minimum  error  of  probability  because  they  fail 
to  consider  the  nonlinear  phase  terms. 

Nielsen  integrated  numerically  with  a constant  one  sec- 
ond interval  in  increasinp  increments  alonp  the  pulse  re- 
sponse until  a time  delay  was  found  that  yielded  a maximum 
output.  However,  the  time  delays  that  Nielsen  obtained  for 
his  filter  are  suspect  by  observation  of  his  Fipure  4 (Ref 
3i20).  It  is  surmised  that  Nielsen  had  an  error  in  his 


profTam  which  resulted  in  a narrow  v/indow  of  integration 
time  rather  than  the  one  second  interval.  This  is  because 
his  delay  times  correspond  to  the  peak  of  the  pulse-response. 

By  duplicating  the  method  tliat  Nielsen  vised,  corrected 
values  of  time  delays  for  this  filter  wore  obtained.  The 
results  are  shown  in  Table  I. 


Table  I 
Filter  Delay 


BT 

T, 

N ielsens 

sec 

Corrected 

2.0 

* 

0.38 

1.5 

0.65 

0.52 

1.0 

1.0 

0.74 

0.75 

1.04 

0.5 

2.0 

1 .60 

♦Not  computed  (From  Ref  3«i^) 


BPSK 

The  probability  of  error  equation  for  BPSK  with  no  in- 
tersymbol interference  is  first  developed.  .M though  this 
result  is  well  known,  the  development  serves  as  a blueprint 
for  BPSK,  QCM  and  QP3K  in  the  presence  of  intersymbol  inter- 
ference . 

Transmitted  BPSK  sin'‘als  may  be  represented  as 

x(t)  = Vp  2 m.(t-iT)  cos  ('.it  + 0(t))  (3) 

i ^ 
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where  m^(t)  is  a T-second  rcctan^’iular  pulse  of  unit  ampli- 
tude with  si^n  - 1 deiiendir.;-;  on  the  i^^  data  value,  and  9(t) 
is  a random  process  representinr  phase  instabilities.  For 
this  th"  's,  T will  always  be  one  second.  The  received  si^'- 
nal  is  y(t)  = x(t)  + n(t),  where  n(t)  is  zero-mean,  white- 
Gaussian  noise  with  power  spectral  density  No/2  (watt/Hz) 

(Ref  4i  Chapters  5t6)  or  (Ref  15«  Chapters  1,4). 

The  baseband  product  of  y(t)  and  the  cosine  reference 
signal  r(t)  after  lowpass  filtering:  is 

y’(t)  = \f?  2 m(t-it)  cos'}(t)  + n (t)  (4) 

i ^ 

where  ^>(t)  is  the  phase  error  between  the  incoming:  phase 
0(t)  and  the  estimate  of  0(t)  made  by  the  VCO. 

The  i^^  data  siy;nal  out  of  the  inte.'trate  and  dump  (l/D) 
circuit  is  d^  = (Odt,  where  I is  the  data  time  iT  + t 

to  iT  + T + T seconds.  For  this  thesis,  the  data  will  be 
obtained  in  the  0^^  data  interval.  Assuminp  that  the  0^^ 
data  bit  m^  = 1 , and  that  4(t)  is  constant  over  the  inte^Tci- 
tion  interval,  the  output  from  the  I/D  circuit  is  = 

T\fp  cos(b  + N,  where  K is  a Gaussian  random  variable  with 
zero  mean  and  variance  NoT/2 . 

The  probability  of  a decision  error,  based  on  the  hy- 
pothesis that  the  0^^  data  bit  is  one  and  the  phase  is 
constant,  is  P(error|mQ  = 1,®)  = ?(T\^  cos4»  + N < 0).  But 
this  is  just  the  "tail”  of  the  Gaussian  density  of  the  ran- 
dom variable  N integrated  from  negative  infinity  to  -TV^  cos$. 

In  terms  of  an  error  function  complement  I 
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P(error|mp.  = orfc(\|—  cooij)  ( 

^ V No 

The  total  probability  of  error  conditioned  on  the  phase 
error  is  then 

P(error|$)  = - erfc(^j;j^  cos<I>)  (6 

for  equal  a priori  probabilities. 

BPSK  with  Intersytnbol  Interference 

In  the  previous  section,  the  cotiditional  probability 
of  error  equation  was  developed  for  BPSK  whore  the  trans- 
mission filter  was  assumed  wide  relative  to  the  data  band- 
width. This  assumption  is  no  longer  valid  as  DT  products 
become  smaller.  The  modulation  pulse  term,  m^(t-iT), 

is  distorted  and  spread  in  time  as  it  passes  through  the 
filter. 

The  input  to  the  data  detector  with  m interferin.n  pulse 
is 


y'(t)  = Vr  Z p(t-mT)  cos(;ot  + 0(t))  + n(t)  (7 
m 

where  p(t)  is  the  pulse-response  of  the  filter  convolved 
with  the  0^^  data  pulse.  This  expression  shows  the  addi- 
tive nature  of  the  intersymbol  interference  relative  to  the 
desired  data  epoch  (m  = 0). 

With  the  ni,pn  of  p(t)  held  positive,  correspond  in/'-  to 
a one  beinp:  sent,  the  noiseless  output  of  the  I /D  circuit  i 

S = VFT  cosCC  X |p(t)|dt  + Z X^p(t-mT)dt7  (B 
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where  the  last  term  represents  an  infinite  sninmation  of  the 
interferers  seen  in  the  desired  data  epoch,  and  I is  the  0^^ 
interval  from  t to  r + T . 

For  computational  purposes,  a trvincation  of  the  m in- 
terferinf;  data  pulses  is  necessary.  The  number  of  future 
pulses  (or  bits)  are  determined  by  the  ^oup  delay  of  the 
filter  for  a desired  BT  product.  Past  bits  are  conti/^uously 
included  for  the  worst  3T  case  until  a power  threshold  rela- 
tive to  the  desired  data  pulse  is  reached  (20  dB  down). 

Table  II  shows  the  past  and  future  bits  considered  in  the 
numerical  computations  for  DFSK  and  QCM.  The  number  of  tits 
for  the  QPSK  case  is  discussed  in  that  section. 


Table  II 

Number  of  Bits  (K3) 


BT 

Past  Bits  (KP) 

Future  bits  (KF) 

0.5 

7 

2 

0.75 

7 

2 

1.0 

1 

1.5 

8 

1 

KB  = KP  + KF 


The  probability  of  error  for  5PSK  conditioned  on  a 
phase  error  and  a truncated  data  sequence  is 


1 \hib 

P(el$,data)  = - orfc(\J — coc.^(S  lp(t)  dt  + 

2 1 No  I ‘ 

KB 

Z X-r  p(t-  m T)dt) ) (9) 

m^’o 
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where  the  truncated  summation  includes  only  those  KB  in- 
terferinf;  bits  from  Table  II. 

With  the  weak  ass'umption  that  intorsymbol  interference 
and  phase  error  are  independent,  the  probability  of  error 
conditioned  on  the  phase  error  is  approximated  as 

2KB 

P(e|$)  = S P(ej<S  ,data.  ) (10) 

2KB  i=i  1 

where  the  conditioned  error  equation  is  summed  and  averafted 
over  all  possible  bit  combinations  for  a data  stream  of 
lenfcth  KB  + 1 . 

QCM  with  Intersymbol  Interference 

The  received  sii-rnal  for  QOM  is  an  altornatinf  auadra- 
ture  sifjial  that  is  binary  modulated.  The  .JCM  signal  is 
best  represented  as  a BPSK  sipnal  in  which  the  carrier  is 
shifted  by  tt/2  radians  every  datum  pulse. 

The  input  to  the  data  detector  is 

y'(t)  = \/p[  Z p(t-mT)  cos(wt  + 0(t))  + 
m 

even 

2 p(t-mT)  sin((i't  +e(t))]+n(t)  (11) 

m 

odd 

where  p(t)  is  the  pulse  response  rcsultinf:  from  the  convo- 
lution of  the  desired  datum  pulse  m(t)  with  the  impulse- 
response  of  the  transmission  filter  as  in  the  previous  sec- 
tion. The  receiver  of  Fi/ture  3.  pare  8,  has  a clocked 
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reference  si^Tial  so  as  to  always  output  the  data  on  the  co- 
sine lef:  of  the  Costas  loop;  exact  timinp  for  the  clock  is 
assumed . 

Aptiin,  with  the  si('7\  of  p(t)  held  positive  oorrespond- 
inf,  to  a one  boinp:  sent,  the  noiseless  output  of  the  l/D 
circuit  is 

S = \fPTlcov.ii[S.  |p(t)|dt  + £ i‘jP(t-mT)dt)  + 

even 
m / 0 

sin4>(£  / pCt-mTldt)"!  (12) 

m ^ 
odd 

This  equation,  when  compared  to  Eq  8,  shows  the  reduction 
of  the  odd  interferers  relative  to  the  desired  datur.  epoch. 
This  result,  when  compared  with  ilq  15  in  Nielsen's  thesis 
differs  in  the  addition  of  the  sinC  torn.  This  te:-m  was 
included  because  it  represents  the  demodulation  distortion 
caused  by  the  sine  channel  when  perfect  coherent  detection 
cannot  bo  maintained.  As  the  carrier  track  inf':  plutse  error 
increases  to  n/4  radians,  the  QCM  sif^nal  from  the  I/D  cir- 
cuit def.enerates  to  that  of  the  BFSK  case.  The  one-shot 
probability  of  error  for  QCM  conditioned  on  a particular 
phase  error  and  data  sequence  is 

P(e|4),data)  = ~ erfc[  { cosi>(/j  ) p(  t)  ] dt  + Z /j,  p(  t-mT)dt) + 

even 
m / 0 

sin4>(Z  /t  p(t-mT)dt)}1  (13) 

m ^ 
odd 
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"SK  with  T !'. t tM’ 


As  stK>w!i  it)  V' l.'Mit'i'  \,  s,  ^i'SK  cat)  be  t'cci'cscatcvl 

as  two  li)iic)'t't)dcn t bittasy  chaaacla.  'I’hc  pfoba'c i 1 i t y o'."  a 
coiTcct  b.t'v' is.  ioi)  is  I'vd^')  - ( 1 - i' i | ^ d v 1 -i' ) | .1  )\  wnese  c. 

and  i\,  arc  the  independent  cosine  atid  sine  channel  errors, 
and  where  the  phas.e  error  is.  the  s..)!ne  for  both  channt'ls.  i’.\et' 
In  wi-itiny,  this  exprt'ssion,  it  is  assumed  th.it 
the  declttion  varia'.'les  t'rom  the  sint'  or  cosine  I ci’scaits 
(d^^  or  d^,  I independe:)t  . Tltis.  is  valid  sit'.ce  t!'.e  nois-e 

variat'Ies.  ar)d  N,  are  nrtcorre  la  ted.  daussi.m  random  vn:i- 
ahles.  j het)Ci' , tl)ey  are  also  independent. 

The  sis.nal  input  to  the  data  detector  is. 

y ' ( t ) ■-  yT[S  p , ( t- ml ) cos.  (:i' t + t>  ( t ' ' ■< 


Z p,(t-mlO  sia(,i't  ) t ) ^ •>  ntti  tl';'' 


With  t!)e  sist)  ot’  the  0^'  cosir.e  datnm  p.ilse  hold  posi- 
tive, the  oatput  from  tl)e  cosine  1 'd  cit'cait  is 

- \/r'T[cos.:(j'  [p,(t)|dt  f z / i\,(t-mddt  + 

ml 


sin4>(Z  S it-inT\it''’  + N, 

ml  "’  1 


where  tl)e  first  term  is  the  nois.eless  s.ip,T)al  from  tb.e 
cosine  I I''  circuit.  A I tt'rn.ite  I y , with  li\e  siyn  of  ths'  0^^ 
sine  datvim  pvilse  held  positive,  d,  - S,  + N,,  .iT)d  N,  .vnd  N, 
are  IndeptMuli'n  t ilaiis.sian  nois.e  t.'sms.  w)th  oero-me.i!)  .vnd 


variance  NoT/2. 


Writin,‘"  the  modulation  set  for  the  two  independent 
channels  as  t 1 and  t 1.  where  the  subscripts  c and 
s refer  to  the  cosine  or  sine  channel,  the  probability  of 
an  error  for  the  cosine  channel  is  r(e^|<I»)  = 


P(e|m^  = +1,4,  data)  + - 1 )P  (e  | rr.^-  - 1 ,$  , data ) . VJith 

m and  m both  beinp  identically  distributed,  the  one-shot 
c s 

probability  of  an  error  is 


P(e^(4,data)  = - erfc[||^  Cccs4(w'' 


dt  + Z 
m 

0 


s 


Pc(t 


•m'r)dt)' 


sin4i(Z  J*  p (t-mT)dt)]'!  (l6) 

ra  I s 


Averaf^inr;  over  all  possible  values  of  a truncated  data  se- 
quence for  both  channels. 


1 

7 


2KB+1 

Z p(e^|  4, data.  ) 

1 ^ ^ o'  1 


(17) 


which,  by  symmetry,  is  also  the  probability  of  error  for 
the  sine  channel.  Psinp  the  equation  for  a correct  deci- 
sion, the  probability  of  a total  symbol  error  for  QFSK  con- 
ditioned on  the  phase  error  is  P(el4)  = 2P(e^j 4 ) -F^ (e^l 4 ) . 

In  order  to  make  a fair  compari.scn  between  QFSK  and 
QCM,  the  symbol  error  for  QFSK  must  somehow  be  related  to 
the  s^Tubol  error  for  QCM.  If  the  QFSK  symbols  are  encoded 
with  a Gray  code,  then  the  bit  error  probability  for  QFSK  is 


P(ejjl4)  = F(e|4)/lor2N 


(18) 
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whore  N ii>  nur.hor  of  ^iFSK  I'.ymbolo  (four)  and  the  aynbol  error 
ip  small  (:\ef  Since  the  symbol  error  for  ddif  is 

also  t!;e  bit  error,  a common  detiominator  t^ow  exists  for  a 
comparison  of  the  two  systems. 

In  order  to  nnmorically  compute  Sq  17  in  a reasonable 
amount  of  time,  the  number  of  truncated  bits  coiisidei'ed  in 
a channel  must  be  reduced  from  those  in  tb.e  binary  cases. 


idered 

in  Table 

11 

were  determined 

u.a  t ion 

(PT  - 0. 

o . 

Inspection  of 

the  wor 

s.t  case 

f e r 

dPSK  (PT  = 1.0^  i 

shows  that  ir'. terferers  witliin  10  dP  cati  always  be  iticluded 
for  this  filter  when  one  future  bit  arid  three  past  bits  are 
cons  idered . 

PFSK  Carr  i-'r  Track  in-"  ret~fermar.ce 

For  a received  ITSK  sijTial  plus  noise  as  riven  in  Cvpaa- 
tion  d,  the  inpvits  to  the  multiplier  prior  to  tb.e  loop  fil- 
ter in  Fit-t.  d are 

A(t)  cosi  n^(t)  (cosine  lee:) 

A(t)  sini  n.,(t)  (sine  ler.)  (Id) 

where  nj(t)  and  r..^(t)  are  c.ero-meafi  iridepenaent  daussiart 
noise  process.es  with  vari.ance  No:-.  Ar,.i  .\{l)  ~ '^'T'Zp  ( t-mT  ^ is 
the  sum  of  all  interfer in,-’:  pulses  "ind  data  seen  in  the  de- 
sired integration  time,  attd  is  ;\ssumed  to  pass  throue;h  the 
lowpass  filters  in  the  loop  with  no  distortion.  Tlie  two- 
sided  bandwidth  of  the  lowpass  filter,  OP,  is  .issumed  to  be 
one  Hz  so  as  to  correspond  to  the  banviwidth  of  .v  matched 


2? 


IL 


filter  for  one  eeoor'.d  pulrer. . 

The  or.tput  f’-om  the  r.viltiplier  applied  to  tV.e  loop 
filter  la 

( 1 1 

e(t)  - — L>iadw(t)  + A(t)  coat  ( t , ( t ) + 


A(t)  ainC(t)n,(t)  + r.,(t)n., (t) 
1 1 *• 


(t) 

are 

iridepeadei'. 

t of  each 

o ti'.er 

pro 

aea  t 

V'  ^ ■'1  ^ 

U » i (.  V.'  V 4.1  k 

e f fee  ti ve 

ao  iae 

0 f 

ad.d 

itive  Gauaa 

iar.  t'.oiae 

ar.d 

interayn'eo I iaterferer.ee  applied  to  Die  loop  with  aa  auto- 
correlat iofi  of  (Sef  iiq 


R '(t)  = R + r;;  (t'i 

^ ^ II  J • V 4 * ^ 


For  phaae  erroi'a  tl'.at  are  ar.all  (<d0°',  aii'.fC'  of  the  aie;aai 
term  ia  well  approx  iria  ted  hy  f v , arui  the  Coataa  loop  oaa  he 
reprer.ented  aa  a 1 iitear  ayatem  aa  aiiowr.  ia  Fipaire  , pa.ae  t d- 
Inherer.t  ia  writia,-.  th.e  effeotive  tioiae  tej-m  at  th.e  iriput  of 
the  Coataa  loop  are  aaaamptioaa  ootioerain,''  the  t\oiae  atid 
phaae  error  prooeaaea.  Ta  eenera^l,  the  rioiae  proeoaa 
ia  a functioa  of  t''.e  phiaa-'  e:a-o’a  ^(t):  'aowever,  with  il'(t) 
asavimed  to  he  alowly  varyir,.'  with  time  (i.e.  coaataat  over 
several  data  iatervala',  a^(tl  a.ad  (*(t)  oaa  he  ahowti  to  he 
effectively  iadepeadent  of  eac’e.  other  (Ref  Id  if  St)  aad  (Kef 
h'q  h-dfl'.  Thia  allows  the  to*al  effective  aoiae  a'(t^ 
to  be  super impoaed  hack  through  the  amplifier  with  :'a\  appro- 
priate ritin  at  the  iaput  of  the  loop. 
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Now  A‘'(t)  ropror.onts  the  tir?.o  a.vera;',inf.  of  thr  error 
sifTial  over  the  coherence  tine  of  the  loop  filter.  For  the 
simulaticn  of  the  loop  filter  in  this  thesii?,  the  coherer.co 
time  of  the  loon  filter  was  chosen  to  be  ter.  seco.nds  . This 


corresponds  to  usine  ten  data  pulses  to  contribute  to  the 
error  sifTial  that  drives  the  VCC. 

The  variance  of  the  phase-error  de.e  to  the  intersynbol 
interference  and  additive  haussian  noise  and  conditioned 
on  a particular  data  sequence  is  appro.x i.r.ately 

2 2BtS^»(0) 

C.data  = — " ' _ — (2i 

(A^)" 


where  ci'-'-ed  loop  banuwidt!;  and  the  power  spectral 

density  of  the  r'.oise  Svf'  is  ass’uned  cor.stant  ovei'  (lef 

ij 

I4i251  - 253) • This  assumes  that  the  closed  loop  bandwidth 

is  much  narrower  tha.n  the  bandlimited  .noise  processes 

n,  (t)  arid  ( t)  . 

1 

Computation  of  fhe  time  average  of  the  input 

sifTial  to  the  loop  filter  is 


A'^Ct)  = ^ Z p‘-(t-mT)dt  ^ /‘L  V v p(  t-mT)p  ( t-ivT)  dt^  ( 2 ^ ) 

L 0 m n 

0 m /n 


where  the  received  siftnal  power  F is  assu.med  constant  over 
the  coherence  time  of  the  loop. 


From  this  equation,  it  can  be  scon  that  the  net  effect 
of  ISI  is  to  reduce  the  averase  do  sis.nal  that  is  applied 
to  the  VCO  and  to  cause  a random  variation  of  the  error 


t 
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signal . 


Computation  of  Applying  the  Fourier  trar^sform 
to  Eq  21,  the  approximate  power  spectral  density  of  the  ef- 
fective noise  seen  by  the  loop  is 


=n- 


(0)  = HS 


R^(0) 


(S|)2 


(24) 


where  is  the  total  average  pov/er  or  the  expression  for 

2 

A (t)  that  was  developed  in  the  preceding  section. 

The  phase  error  variance  conditioned  on  a particular 
data  sequence  and  due  to  I3I  and  additive  Gaussian  noise 
is  rewritten  as 


1 No-. 

^4>,data  “ -J 


PI- 


where  I = CJ'(^i'Zp^(t-mT)dt  + 2p(  t-mT)p(  t-nT)dt']  (25) 

O m G m n 


m n 
m=n 


Averaging  this  conditional  variance  over  all  of  the  possible 
combinations  of  data  pulses  in  the  coherence  time  of  the 
loop. 


? 


CTj/Tsl 

.2 

°C>,data; 


(26) 


where  CTj/Tg3  is  the  greatest  integer  value  of  the  effec- 
tive number  of  symbols  used  by  the  loop  in  making  a phase 
estimate.  For  no  ISI,  I is  just  o.ne  and  Eq  25  reduces  to 
Eq  8. 3 •13  as  given  in  Stifflcr  (Ref  \hi2k7). 
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acM  Carrier  Trackinr  Ferformanoe 

For  a received  QCM  sirnal  plus  noise  as  fiiven  in  Hqua- 
tion  11,  the  inputs  to  the  multiplier  prior  to  the  loop 
filter  of  Fipiire  3 

Ae(t)  cos4»(t)  + Ao(t)  nin^>(t)  + n^(t)  (cosine  lep:) 

Ae(t)  sin<5(t)  + Ao(t)  cos4i(t)  + n^{t)  (sine  lop)  (*37) 


whore  Ao(t)  and  Ae(t)  are  the  summations  of  the  odd  and 
even  interforinp  pulses.  The  output  from  the  multiplier 
applied  to  the  loop  filter  is 


e(t)  = (Ao^  ( t ) + Ae'^  ( t ) ) sin2C)(t)  + (Ae(t)  sin(I»(t)  + 

2 

Ao(t)  cosC-C  t)  )n^  ( t)  4 (Ao(t)  3inC>(t)  + 

Ae(t)  cos$(  t)  )n-,(  t)  + n,(t)n.,(t)  (28) 


where  C(t)  and  n.(t)  and  n^vt)  are  independent  as  before. 
Also,  Ao(t)  and  Ae(t)  are  assumed  to  bo  independent  since 
they  are  formed  from  a sei’ies  of  independent  data  pulses. 
The  last  three  terms  represent  the  total  noise  oit-;;nal  n'(t) 
which  cor\sists  of  additive  Caussian  noise  and  intersymboi 
interference  .and  has  .an  .autocorrelation  of 

The  autocorrelations  of  Ao(t)  and  Ae(t)  are  e.asily  ob- 
tained from  the  autocorrelation  of  A(t)  for  the  PFoK  case. 
The  squ.are  of  A(t)  is  written  as  the  svim  of  the  odd  and 
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2 2 2 

even  terms  squared,  c.f.  A (t)  =Ae  (t)  +Ao  (t)  + 

2Ao(t)Ae(t).  Therefore,  the  time  avera^v  of  the  input 

sipTial  to  the  loop  filter,  Ao2(t)  + Ae^(t),  is  just  A‘^(t) 

from  the  BFGK  case  minus  all  of  the  terms  of  odd  times 

even  products  in  the  second  intepral  of  Eq  23.  Likewise, 



the  total  averapie  power  CRAe(^)  RAo(0)1  is  Ao  (t)  + Ae‘'( t ) . 

Trackinr  Penal  ty  for  QPSK  versus  PPoK'/ICm  with  no  1ST 

Althoufth  an  expression  was  not  obtained  for  the  phase 
error  variance  of  an  N = 4 Costas  loop,  the  ent'rfcy  pene.lty 
for  tracking  a quadrature  versus  a binary  phase  shift  sif’:- 
nal  with  no  ISI  can  be  used  as  a desipn  cons ideration . The 
motivation  for  developinf  this  expression  is  as  follows: 
Since  the  error  performance  for  QCM  is  very  close  to  that 
of  QPSK  for  the  same  data  rates  and  transmitter  power,  per- 
haps QCM  could  be  used  to  an  advanta,te  whore  phase  instabil- 
ities are  a problem.  The  bandwidth  of  a QCM  carrier  traok- 
inp  loop  can  be  increased  so  as  to  obtain  a faster  response 
time  to  these  phase  changes  than  can  be  obtained  from  QPSK 
for  a piven  transmitter  power.  Alternately  stated,  when 
the  maximum  phase  error  and  response  time  of  a system  is 
specified,  what  is  the  increase  in  Sb/rJo  that  must  be  ap- 
plied to  a QPSK  system  so  it  has  the  same  track in.y  charac- 
teristics as  a QCM  system? 

Equating:  the  expressions  for  phase  error  variance  for 
BPSK  (also  QCM)  and  QPSK  under  no  ISI  yields 

TsBL[l/(Z-a)  + l/2(Z-a)2l  = T’^P'^Cl/Z  + 9/2  Z"  + 

6/z3  + 3/2  Z^ 1 (30) 
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where  Z is  the  bit  i3 i(''r\al-»merr,y- to-no itv.'  ratio  for  ^loK, 
Tg  and  T'^are  the  r.ymbol  Intorvalr. , and  and  H'^^ai'e  the 
bandwidtlu;  of  the  cloned  loopr.  (Hof  l-’tt  IPi  8.T.13).  For  a 
fair  comparison  of  the  two  syi  toms  the  infoi-ma t ion  rater, 
are  the  aame  (T',  = To/Z),  and  the  responr.er.  of  the  track- 
inf:  loops  are  the  same  (8%  " P..  ) . 

Xj  xj 

Fixinp.  the  ripht  hand  side  of  liq  TO  at  sonu'  si,qavl- 
enerpy- to -no ice  ratio  Zo  and  then  makinp,  the  chanpe  cf 
variable  Z - a - X produces  a quadratic  e:)uation  ir.  X. 


A(Zo)  X"  - X 


2 = 


(31) 


where  A(Zo)  is  tne  ripht  hand  side  of  L'q  ia  brackets 
evaluated  at  Zo.  Pe cause  the  enerpy  sipnals  cannot  ne 
nepative  and  the  discriminant  is  always  preate:-  than  -/.ero. 
the  sipnal- ener;-y- to- no ise  ratio  of  ti'e  QCM  case  ir, 


1 + Vl  -»  2A(Zo’) 


X = 


for  the  same  phase  error  variance,  response  to  phase  insta- 
bilities, and  information  rate  as  the  QPSK  ease. 


III.  Numerical  I^caul ts 


In  Chapter  II,  equations  are  riven  for  the  conditional 
probabilities  of  error  for  BFSK , wCM,  and  QPSK  in  the  pres- 
ence of  intersymbol  interference.  Also  riven  are  the  aver- 
ape  phase  error  variances  for  BPSK  and  QCM,  and  the  onerry 
penalty  for  trackinr  a quadrature  versus  a binary  modulated 
sirnal.  In  this  chapter,  numerical  results  from  those  equa- 
tions are  obtained  in  order  to  compare  QCM,  BPSK,  and  QPSK 
systems.  First,  the  philosophy  for  includin.r  only  those  Iri- 
terferinr,  pulses  within  a certain  power  threshold  is  discus- 
sed, and  a comparison  with  Jones's  results  is  shown.  The 
conditional  error  probabilities  are  thoii  compared  for  PP5K, 
QCM  and  QPSK.  Following:  that,  the  phase  error  variances 
for  BPSK  and  QC’.l  are  examined.  Finally  the  enor^-:y  penalty 
for  track  in,";  a QPSK  si^tnal  versus  a QCM  or  BPSK  sipnal  is 
presented.  The  overall  conclusion  is  that  QCM  performs 
better  than  BPSK  and  is  only  sli,"htly  inferior  in  data  de- 
tection performance  to  QPSK. 

Truncation  of  Data  Sequence 

A major  problem  in  this  thesis  was  how  many  interferinp 
pulses  should  be  included  to  obtain  sufficiently  accurate 
numerical  results.  Jones  has  stated  that,  except  for  low 
BT  products,  satisfactory  results  can  be  obtained  in  eval- 
uatinp;  error  probabilities  by  considerin,"  only  the  most  ad- 
jacent interferer  from  the  past  data  bit  (Ref  llil26).  How- 
ever, ho  did  not  state  at  what  BT  product  this  assumption 


becomes  invalid.  The  number  of  past  bits  (KP)  used  for  the 
binary  modulation  schemes  for  this  thesis  was  chosen  to  be 
KP  = 9 - KF,  where  the  future  bits  (KF)  are  determined  from 
the  {’iroup  delay  in  Table  I,  Chapter  II.  V.'ith  this  ali'orithm, 
the  power  levels  of  all  pulses  outside  of  the  truncated  data 
sequence  are  at  least  20  dB  below  the  power  of  the  current 
datum  pulse  for  a PT  of  0.5  (the  worst  case  of  ISI).  While 
this  results  in  including  terms  below  the  20  dP  threshold 
for  hifher  PT  products,  their  inclusion  should  only  serve 
to  increase  the  accuracy  of  those  results. 

The  power  threshold  of  20  d3  was  verified  by  comparing' 
the  conditional  probability  of  error  for  BFSK  for  two  ad- 
jacent interfering’;  pulses  and  for  the  nii;e  ii^terfer i.r,-’;  pulses 
from  Table  II.  These  results  are  shown  in  Table  III  for  a 
bit  s ipnal- cner,ty -to -no ise  ratio  of  10  dP.  This  comparison 
shows  the  necessity  for  includin,t  more  than  two  adjacent 
bits  as  the  PT  product  approaches  O.S  to  obtain  more  than 
one  significant  fipure  of  accuracy.  The  -20  dB  power  thres- 
hold is  assumed  to  be  adequate  with  the  following  reasoriir',f: « 
The  two  adjacent  bits  for  a PT  = 1.0  correspond  to  includin;'; 
all  pulses  within  a power  threshold  of  only  -15.5  dP.  In- 
creasing the  number  of  bits  to  nine  for  a PT  - 1.0  results 
in  considering  all  pulses  within  a threshold  of  -42  dPj 
however,  this  only  changes  the  fourtl\  significant  dipit  in 
the  probability  of  error  results.  Thus,  i power  tiireshold 
of  -20  dB  for  the  worst  case  ( PT  - 0.5)  "'as  arbitrarily 
chosen  to  lie  between  those  vaUies,  and  is  expected  to 
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Table  III 


BFSK  Conditional  Probability  of  2rror 
(Fliaae  Error  = 0°) 


BT 

Two  Adjacent 

Fulses 

KP  ^ KF  = 1 

P(el*) 

N ine 

Fulses 

KP  + KF  ^ d 

KF 

KP 

0.5 

1.6730(-3) 

2.0701 (-3) 

L, 

7 

1.0 

5.1155(-5) 

5.1l3l(-5) 

1 

6 

1.5 

2. 1854 (-5) 

2.1d43(-5) 

1 

8 

Note*  The  nof'n.tivo  ir.te,-''er  in  paronthenon 
follow iiip;  oL-ch  entry  repi-  en'^ntr  the  powi'*r  ot 
ton  by  which  the  enti'y  should  bo  ;;iul  tipi  ied . 


produce  probability  of  error  results  that  are  accurate  in 
the  third  significant  di,tit.  Similar  reasoninf  was  used 
in  the  QFSK  case;  however,  the  threshold  was  arbitrarily 
established  at  -30  dF. 

The  BFSK  cane  for  no  phase  error  is  compared  to  Jor.es 's 
Fif;.  6.  This  comparison  is  shown  in  Table  IV  as  an  increase 
in  the  sipnal  enerry  (62b)  from  the  ideal  case  of  no  ISl 
(BT  = <*')  required  to  maintain  a bit  probability  of  error  of 
10"^  with  perfect  coherent  detection  assumed.  It  sho\jld  be 
noted  that  Jones  defines  his  FT  on  the  basis  of  a two-sided 
bandwidth;  therefore,  his  FT  of  2.0  corresponds  to  a FT  of 
1.0  as  used  in  this  thesis.  This  comparison  was  important 
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for  three  reasons.  First,  the  correlation  of  the  two  methods 
served  to  verify  that  the  pror,ra;ns  in  this  thesis  were  oper- 
atinf:  correctly.  Second,  the  increas  in.-;  difference  be  tweets 
the  two  methods  as  the  ?T  is  decreased  verified  Jones's  con- 
clusion that  more  than  one  adjacent  bit  is  required  for  low 
BT  products.  Finally,  for  3T  levels  of  1.0  or  less,  more 
than  one  adjacent  bit  should  be  included  to  obtain  suffi- 
ciently accurate  results. 


Conditional  Frebabil  ity  of  ET-ror 

The  conditional  error  probabilities  of  PFSK,  aiid 

QFSK  avera^'ed  over  the  truncated  data  stream  and  cot\d  i t iot\ed 
on  phase  errors  of  0*^,  10*^,  and  20°  a.re  compared.  Thesi 
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conditional  probabilities  are  also  examined  relative  to  the 
ideal  signal  case  (no  in ters^mibol  interference  and  perfect 
coherent  detection). 

BFSK  versus  QC'I.  In  Fifuros  7,  5,  and  d the  conditional 
error  probabilities  of  the  BPSK  and  QCM  sipjials  with  phase 
errors  of  0°,  10°,  and  20°  are  shown.  These  probability  of 
error  plots  show  the  improvement  over  normal  BPSK  by  usinp 
QCM.  The  defxadin^  effects  of  an  increasin.t  phase  error  are 
apparent  when  those  figures  are  compared  to  each  other. 
Tables  V,  VI,  and  VII  compare  BPSK  and  QC?.:  for  selected 
sif:nal-enerf:y-to-noiso  ratios.  The  tables  show  that  al- 
though QCM  has  a lower  probability  of  error  than  BPSK,  it 
is  more  sensitive  to  a chanf:e  in  phase  error  than  BPSK. 

Nielson  compared  these  two  modulation  techniques  v.'ith 
an  improvement  factor  defined  as  the  Eb/No  (dB)  a QCI-i  sig- 
nal requires  to  obtain  an  error  probability  of  lO"^  minus 
the  Eb/No  (dB)  a BPSK  sif:nal  needs  to  produce  the  same  prob- 
ability of  error  (Ref  3*29.  Table  VI).  He  found  an  extrem- 
ly  hif:h  improvement  factor  (jxoater  than  26  dB)  with  a BT 
of  0.5.  This  author  feels  these  are  a result  cf  the  in- 
correct pq'oup  delay  and  tne  missintt  term  in  Nielsen's  Eq  16. 

For  this  thesis,  an  improvement  factor,  3.  is  defined 
as  the  difference  in  Eb/No  (dB)  for  two  sirnallinf;  schemes 
in  order  to  produce  the  same  specified  probability  of  error. 
A power  penalty  is  defined  as  the  additional  power  (dB)  be- 
yond the  ideal  case  (BT  - '»)  required  to  maintain  the  sane 
specified  probability  of  error.  In  Table  VIII,  this  im- 
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Fir^  BPSK  ami  ^Cy.  Bit  Frror  r.Tl'orT.anoo  (Phase 


= 10°) 


nrror 


Fif-'.  9«  PF3K  and  QCM  Bit  i£rror  Forforr.ance  (Phase 
Error  = 20°) 
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Table  V 


Conditional  Probability  of  Error 
For  a FT  of  0.5 


Eb/l^o 

I 

1 

Phase  Error, 

F(eU) 

dB 

de  frees 

BPSK 

3CM 

0 

0 

0.1375 

! 0.1280 

0 

10 

0.1400 

0.131^ 

0 

20 

0.1514 

0.1440 

6 

0 

0.236o(-i)  ; 

i 0.1l3l(-l) 

6 

10 

0.2S05(-1) 

1 O.n2o(-1) 

6 

20 

0.2P6F(-1 ) 

0.l83-i(-l) 

10 

0 

0.2070(-2) 

0.  lo0'3(-‘!) 

10 

10 

0.2316(-2) 

0.2?00(-'^) 

10 

20 

0.3210( -2} 

C. 6634 (-3) 

12.5 

0 

0.12Sl(-3) 

0.11  S6(-S) 

12.5 

10 

0.1532(-3) 

o.2o5:(-s) 

12.5 

20 

' 0.2535(-3) 

. 1 

o.nr9(-4) 

Note  I The  nerativ'e  intewr  in  parenthot^is  follow- 
ing: each  entry  in  the  table  rerrooenta  tb.e  power 
of  ten  by  which  the  entry  should  be  multiplied. 


Table  VI 


Conditional  Probability  of  Error 
For  a 3T  of  1.0 


b/No 

d3 

Phase  Error 

degrees 

0 

0 

0 

10 

0 

20 

6 

0 

6 

10 

6 

20 

10 

0 

10 

10 

10 

20 

12.5 

0 

12.5 

10 

12.5 

20 

P(el^) 

3PSK  QCM 


Table  VII 

Conditional  Probability  of  Error 
For  a DT  of  1.5 


Eb/No 

1 

Phaac  Error 

P(e|<l  ) 

dB 

defTeoG 

3PSK 

QCM 

0 

0 

0.0239(-1 ) 

0.0!63(-1 ) 

0 

10 

0.9S74(-1  ) 

0.9502(-l) 

0 

20 

0.1062 

0.1056 

6 

0 

0.4320 

0.3800 (-2) 

6 

10 

0.4340(-2) 

0.4307(-2) 

6 

20 

0.6731 (-2) 

0.62o4(-2) 

10 

0 

0.2194 (-4) 

0.1 200 (-4) 

10 

10 

0.2830(-5 ) 

0.1 740 (-4) 

10 

20 

0.591 4 (-4) 

0.4102 (-4 ) 

12.5 

0 

0.30Ql(-7) 

0.8234 (-3) 

12.5 

10 

0.4791 (-7) 

0.l440(-7) 

1 ^ C 

20 

0.169^^(-6) 

0.6803(-7) 

I 
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Table  IX 


Conditional  Probability  of  Error 
For  a QCM/qPSK  BT  of  0. 5/1.0 


Eb/No 

Phase  Error 

P(e|  C*) 

dB 

de grees 

QPSK 

QCM 

0 

0 

0.9571(-1) 

C.12S0 

0 

10 

0.9878(-1) 

0.1319 

0 

20 

0.1083 

0.1440 

6 

0 

0.59l3(-2) 

^ O.llSl(-l) 

6 

10 

0.6563(-2) 

0.1326(-1) 

6 

20 

0.8866(-2) 

0.1834(-1) 

10 

0 

0.5118(-^) 

0.19C7(-3) 

10 

10 

0.64l8(-4) 

1 0.2700(-3) 

10 

20 

0.1237(-3) 

0.6634(-3) 

12.5 

0 

0.1571(-6) 

0.1156(-5) 

12.5 

10 

0.2305(-6) 

0.2457(-5) 

12.5 

20 

0.6993(-6) 

0.1379(-4) 

Table  X 


Eb/No 


Conditional  Probability  of  Error 
For  a QCrVQFSK  BT  of  1. 0/2.0 


Phase  Error 


degrees 


P(e| $) 


0.84S6(-1) 

0.8789(-l) 

0.9731(-1) 

0.3631(-2) 
0.4093 (-2) 
0.5787(-2) 

0.1283(-i|) 

0.i687(-4) 

0.3732(-4) 

0.1065(-7) 

0.1709(-7) 

0.6711(-7) 


0.99o0(-l) 

0.1C30 

0.1139 

0.5UQ(-2) 

C.5751(-2) 

0.8044(-2) 

0.24‘^6(-4) 

0.3280(-4) 

0.7473(-4) 

0.2566(-7) 

0.4426(-7) 

0.2033(-6) 
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Table  XI 


Conditional  Bit  Error  Performance  For 
QCM/QFSK  with  a BT  of  0.75/1-5 
(Phase  Error  = 0°) 


P(el«>) 

EbA'c 

i.dB 

3 .dB 

Power  Penalty 

QPSK 

QCM 

! 

QPSK.dB 

O.l(-l)  1 

4 . ^0 

>1  nn 
^ • ( ( 

-0.16 

0 . ‘■>7* 

0.l(-2) 

7 . .t  Q 

7.75 

-0.26 

O.Q^* 

0.l(-3) 

0 . 1 S I 

0. 37 

-0.21 

0.70* 

0.1(-5) 

11.30 

. i 

11 .46 

- 0 . 1 0 

- 

*Idcal  bit  probability  of  error  for  ^PSK  was  ob- 
tained from  Lindsev  and  Sim.on  (Ref  lOiCir. , Table 

5-3)- 
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rier  trackinr  inatabilitios  are  a problem.  As  an  example, 
in  a low  sifrial- to-noise  environment  where  transreceiver 
instabili oies , Doppler  phase  shifts,  and  random  variations 
in  the  modulating  sipnal  cause  inherent  difficulties  in 
tracking:  the  carrier  of  a quadrature  modulated  sip:nal,  the 
carrier  of  a binary  modulated  sif^nal  can  still  be  tracked. 
This  will  be  discussed  further  when  the  Sb/No  penalty  for 
tracking  a QPSK  signal  versus  a QCM  signal  is  examined. 

Phase  Error  Variance 

The  average  phase  error  variance  is  obtained  by  a nu- 
merical evaluation  of  Eq  2k  for  BFSK  and  QCM  in  the  presence 
of  ISI.  The  phase  error  variance  is  plotted  as  a function 
of  the  signal-energy-to-noisG  ratio  and  filter  3T  products 
times  a closed-loop  bandwidth  in  Figures  14  and  15.  These 
plots  also  show  the  phase  error  variance  for  no  ISI  (BT  = ») . 
For  data  rates  corresponding  to  BT  products  of  1.0  or  higher, 
the  phase  error  variance  is  approximately  the  same  as  that 
for  no  ISI.  For  a 5T  of  0.5  with  sif-qial-energy-to-noise 
ratios  of  6 dB  or  more,  the  phase  error  variance  is  in- 
creased by  at  most  0.18  Rad/'lls  over  the  ideal  case  of  no 
ISI.  The  6 d3  signal  to  noise  ratio  was  a’'bitrarily  chosen 
to  correspond  to  a probability  of  error  of  O.05  for  BP3K 
operating  with  a phase  error  of  CO'^.  This  is  expected  to 
be  the  worst  data  performance  that  could  be  tolerated  with 
direct  data  detection  (no  encoding) . 
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arSK  Tracking  Penalty 

The  bit  si^^ial  enorf;},-  penalty  for  tracking:  a quadrature 
modulated  ^i.qnal  versus  a binary  signal  with  no  I3I  and  the 
same  phase  error  variance  and  loop  response  is  illustrated 
in  Figure  l6.  The  enerfy  penalty  is  defined  as  the  addi- 
tional bit  signal-energjf-to-noise  ratio  Eb/No  in  dB  that  a 
QPSK  system  would  require  as  a function  of  the  3PSK  Et/No . 
While  this  energy  penalty  was  obtained  under  no  I3I  condi- 
tions, it  is  assumed  that  it  can  be  used  as  a rough  approxi- 
mation for  the  tracking  penalty  QPSK  would  suffer  over  QCM 
under  I3I  conditions. 

In  order  to  best  understand  the  importance  of  this  energy 
penalty  consider  the  followingj  In  the  design  of  any  car- 
rier tracking  loop,  the  closed-loop  bandwidth  is  made  as 
narrow  as  possible  so  as  to  mi.nimise  the  effects  of  the 
noise  process  on  the  phase  error,  but  still  wide  enough  so 
the  loop  will  respond  to  the  various  phase  instabilities 
(Ref  15»133"13^’*  In  this  design  process,  it  is  tacitly  as- 
sumed that  the  carrier  tracking  portion  of  the  receiver  has 
a sufficient  front-end  signal- to-noise  ratio  to  maintain  the 
loop  in  a locked  condition.  For  a decreased  input  si^tnal  to 
noise  ratio,  the  closed-loop  bandwidth  cari  be  decreased,  but 
only  at  the  expense  of  reducin.g  the  response  of  the  filter 
to  the  phase  variations.  In  general,  linear  loop  analysis 
predicts  that  this  trade-off  between  the  noise  and  phase 
processes  becomes  more  critical  as  the  nu.mber  of  signalling 
phases  increase  (Ref  lUi2^7)*  In  some  situations,  a low 
received  signal- to-noise  ratio  (SNR)  combined  with  excessive 
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random  phase  processes  will  prevent  proper  loop  track  in,";. 
These  random  phase  processes  are  the  transrece iw-r  instabil- 
ities, Doppler  phase  shifts,  and  variations  in  the  phase 
modulation  (Ref  l')«  Eq  4-8). 

When  improper  loop  trackin^T  occurs  for  QPSK,  ;in  altern- 
ative sipnallin,"  method  such  as  BFSK  can  sometimes  still  bo 
tracked.  As  shown  in  Fi,"uro  t6,  a PPSK  si,";nal  with  an  Eb/No 
of  5-D  dp  corresponds  to  a QPSK  si;7ial  with  an  Eb/Ko  of  7-^ 
dB  for  the  same  information  rate  and  loop  tracking  charac- 
teristics. Suppose,  for  example,  a QPSK  system  wa.s  pla.euetl 
by  phase  instabilities  such  that  it  fell  out  of  lock  with 
an  Eb/No  of  7*  3 dP.  A binai'y  modulated  system  could  bo 
switched  in  at  this  point  and  still  continue  to  operate  ur.- 
til  nearly  5.0  dP  before  it  also  fell  out  of  lock.  Phis 
procedure,  i.e.  switchin,"  in  BPSK  as  a back-up  to  QPSK,  has 
been  proposed  for  DSCS-III.  Naturally,  becausi'  of  the  lower 
Eb/No,  the  probcibility  of  error  will  increase. 

Because  QCM  has  a superior  performance  to  PPSK  and  be- 
cause QCM  is  also  a binary  modulated  si/T.al  it  should  be 
considered  as  an  alternative  back-up  ai,enallin,"  scheme  for 
QPSK.  Clearly,  only  a sli,"ht  improvement  could  be  ,'’ained 
by  usinr  QCM  versus  PPSK  as  a back-up  scheme  where  PT  prod- 
ucts of  1.0  or  hi, "her  are  used.  However,  for  data  rates 
corresponding,  to  a BT  = 0.5  up  to  a 2.7  dB  ,"ain  in  data 
performance  could  be  realised  ever  the  PPSK  system.  Ob- 
viously, no  stron"  conclusions  can  be  made  concernin,"  QCM 
as  a back-up  system  for  QPSK  until  tlio  track  in."  penalty  for 
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Conclusiona 

QCM  faves  better  performance  than  BFSK  when  used  in  a 
communication  system  limited  by  intersymbol  interference. 

For  data  rates  correspondinr  to  BT  products  of  1.0  or  more, 
the  actual  fain  over  BFSK  is  neflifible.  However,  for  an 
increased  data  rate  correspondinc:  to  a DT  of  0.5i  an  im- 
provement of  2.69  dD  is  obtained  for  cin  error  probabilixy 
of  10”^.  Thus,  a chanfe  from  BFSK  to  QCM  could  allow  one 
to  maintain  the  same  probability  of  error  at  almost  double 
the  data  rate . 

QCM  was  found  to  be  inferior  to  QFS’C  when  compared  or. 
the  basis  of  equal  bit  information  rates.  However,  the  dif- 
ference between  QCM  and  QPSK  is  small.  For  example,  QCM 
suffered  at  most  a 0.8  dB  penalty  when  compared  to  QPSK  for 
a BT  pair  of  0.5/1  *0  and  an  error  probability  of  10"-^.  And 
for  a BT  pair  of  0.75/1*5  the  performance  of  QCM  was  almost 
equal  to  that  of  QF3K.  This  almost  equal  performance  suf- 
gests  that  QCM,  and  not  BFSK,  should  be  considered  as  an 
alternate  siftnallinf;  scheme  if  phase  tracking  cannot  be 
maintained  for  QFSK. 

QCM  introduces  no  problems  in  the  carrier  tracking 
capability  of  a Costas  loop,  and,  in  fact,  th'^  average  phase 
error  variance  of  a Costas  loop  was  found  to  be  essentially 
equal  for  BFSK  and  QCM  for  the  BT  products  considered  in 
this  thesis.  For  data  rates  corrospondinf,  to  BT  products 
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of  1.0  or  more,  the  phase  error  variance  is  close  to  the 
ideal  wideband  case  (BT  «) . For  a 5T  of  O.f'i  the  phase 
error  variance  is  increased  from  the  ideal  case  by  at  most 
0.18  rad/Hs  for  received  sirnal-ener^'y-to-noise  ratios  of 
6 dB  or  more . 

Recommendat  ions 

Other  aspects  of  QCM  versus  BFSK  and  QF3K  should  be  in- 
vestigated. A comparison  of  QFSK,  BP3K , and  QCy.  should 
be  made  assuminr  that  the  timing  information  is  in  error. 
This  timing:  error  could  possioly  be  modeled  as  a random 
variation  applied  to  the  eroup  delay,  t. 

Another  area  of  investifntion  is  a comparison  of  the 
three  modulation  schemes  with  an  additional  method  of  ro- 
ducinf:  intersjy'mbol  interference.  As  an  example,  does  the 
application  of  equilization  to  the  three  sipiallinr  schemes 
cause  any  reversal  of  the  hierarchy  of  performance  that  wo.s 
obtained  in  this  thesis? 

A third  area  of  consideration  is  a comparison  of  the 
phase  tracking":  characteristics  of  QCM  versus  QFSK  in  tb.e 
presence  of  intersymbol  interference.  Although  this  thesis 
has  shown  the  Sb/No  penalty  that  a fourth  order  loop  suffers 
over  the  binary  modulated  case  under  no  131,  it  is  not  known 
if  the  effects  of  ISI  would  chanpie  these  results. 
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Appendix  A 


Computer  Frot^rans 

Early  in  the  dovelopricnt  of  the  equations  in  Chapter  II, 
many  common  terms  were  found  to  exist  from  one  equation  to 
another.  For  this  reason,  the  major  program  tasks  wore 
written  as  modules  that  could  be  utilized  by  different  driv- 
er profa'a^s . All  pro^’:rammin/3  was  written  in  Fortran  IV  for 
use  on  the  CDC-6600  computer  system  at  the  Air  Force  Insti- 
tute of  TechnolO£:y,  V.'ri,‘’:ht-Fatterscn  Air  Force  Base,  Chio. 
The  followinf:  descriptions  of  the  modules  or  subroutines  to- 
gether with  the  listings  in  Tables  XII,  XIII,  and  XIV  are 
provided  should  anyone  wish  to  use  those  programs  for  future 
work . 

Subroutine  PULSE 

Subroutine  PULSE  computes  the  pulse-responso  of  an  i^^ 
datum  signal.  Ten  data  intervals  of  the  pulse  response  are 
stored,  at  0.02  second  intervals,  in  the  array  F.  The  sub- 
routine is  used  by  calling  it  with  an  integer  value  (N) 
which  corresponds  to  the  desired  datvim  signal  plus  one. 

Function  SR  computes  the  step-response  of  the  Chebyshev 
filter  in  subroutine  PULSE.  The  function  is  used  by  pas- 
sing it  a real  ar,'niment  which  specifies  the  time  for  the 
quantized  step-response  to  bo  computed  I the  time  values  in 
P range  from  0.0  to  10.0  seconds  in  0.02  second  intci'vals'. 
The  step-response  is  normalized  in  terms  of  BT  products  as 
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discussed  in  Chapter  II. 


Subroutine  DCAV3 

This  subroutine  nvj’iorically  intefratos  any  point  values 
that  are  in  array  PR  and  stores  the  result  in  the  two- 
dimensional  array  FM.  DCAVG  is  used  by  first  declarinf: 
lower  and  upper  limits  for  the  intefration  (BOT  and  TOP). 
These  values  are  passed  to  DCAVG  via  common  block  FOUR. 

The  two  callinr  ar,r;umcnts  for  DCAVG  (I,J)  establish  the 
stora^ie  location  for  the  value  of  the  integral  in  PM. 

Function  FIT . FIT  is  a polynominial  interpolation  rou- 
tine used  in  DCAVG  whenever  a numerical  value  is  required 
by  DCADRE  that  lies  between  the  quantised  data  values.  This 
function  uses  INTERP,  an  Aitken's  Kth  Decree  Polynomial  In- 
terpolation routine  in  the  CDC-66OO  user  library.  The  de- 
piree  of  the  polynomial  used  for  this  thesis  was  three.  This 
choice  was  arbitrarily  made  based  on  the  overall  smoothness 
of  the  pulse-response  in  Fifuro  6. 

Function  DCADRE . DCADRE  pez'forms  the  actual  intofa'ation 
in  subroutine  DC.AVG.  It  is  an  adaptive  Romberf,  extrapola- 
tion routine  of  the  International  Mathematical  and  Statisti- 
cal Library. 

Subroutine  BP 

This  subroutine  computes  all  of  the  combinatorial  possi- 
bilities of  a binary  bit  strin-'.  It  is  used  by  first  com- 
putinr  the  number  of  possibilities  (KNT  - 2-»*M)  where  M is 
t.ne  number  of  bits.  The  subroutine  is  called  with  M and  an 
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integer  from  1 'to  OT  in  steps  of  one.  Therefore,  this 
subroutine  is  called  KNT  times.  The  combinatorial  possi- 
bilities are  computed  as  binary  additions  of  the  previous 
bit  pattern  plus  one.  The  output  is  stored  in  both  IZEP.O 
(corresponding:  to  lof:ical  values  of  '1'  or  '0')  and  IBIT 
(correspondirif:  to  sif:nal  values  of  '1*  or  '-I')- 

Additional  Frorrams 

The  analytical  expression  for  the  stop-respo.nse , as 
discussed  in  Chapter  II,  was  obtained  by  a numerical  analy- 
sis through  the  use  of  FCLY  and  FARTL.  These  are  local  Air 
Force  Institute  of  Technolo^^y  user  programs. 

The  plots  were  obtained  by  usirif:  local  utility  plotting, 
routines  on  a CALGOr'IF  plotter. 

The  error  function  complement  was  obtained  from  SRF . 
This  is  a McLaurin  series  or  continued  fraction  expansion 
routine  from  the  CUC-6600  users  library. 


Prorram  Nota  tlor. 

The  followinr:  is  a list  of  the  major  prorram  variables 
and  arrays  and  their  usaro i 

BT  = bandwidth  tine  product 

TWFI  = two  times  Fi 

C(4,4)  = array  holdinr:  the  coefficients  of  the  Chebyshev 

filter 


X(502) 

Y(502) 

PR (500) 


= arrays  for  the  plots 

= input  array  for  DCAVG  (contains  portions  of  the 
pulse-response  or  products  of  pulse-responses 
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p(500)  = 

PM(10,10)= 

IBIT(M)  = 

TOP 

BOT 

PHE 

TAU 


as  required);  essentially  a scratch  pad  array 
output  array  for  PULSE  (contains  pulse-response) 
output  array  for  DCAVG  (contains  inte^^ral  values) 
bit  patterns 

upper  limit  of  the  integral  in  DCAVG 
lower  limit  of  the  integral  in  DCAVG 
phase  error  in  radians 
delay  times 


Sample  Fro rram , IPSK  Probability  of  Error 

The  program  in  Table  XIII  computes  the  probability  of 
error  for  QPSK  from  Equation  17  in  Chapter  II.  This  program 
was  used  to  obtain  the  plots  in  Figures  10,  11,  and  12  in 
Chapter  III.  The  program  is  very  straightforward.  After 
initialization,  the  outermost  loop  (DO  LOOP  1)  provides 
probability  of  error  curves  for  the  three  BT  products.  First, 
all  of  the  sine  and  cosine  interference  terms,  as  well  as 
the  0^^  datum  pulse,  are  computed  and  stored  in  the  two- 
dimensional  array  FK . Then  nested  DO  LOOPs  (45  and  25)  are 
used  to  compute  all  of  the  combinatorial  possibilities  for 
50  sigr\al-energy-to~ noise  ratio  points;  the  resulting  prob- 
ability of  error  values  for  these  points  are  stored  in  ar- 
ray Y . 


Sample  Program,  3PSK  Phase  Error  Variance 

The  program  in  Table  XIV  computes  the  phase  error  for 
BPSK  from  Equation  26  in  Chapter  II.  This  program  was  used 
to  obtain  Figure  l4  in  Chapter  III.  First,  all  of  the 
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intefrral  values  from  the  I term  of  Equation  25  are  computed 
and  stored  in  the  two-dimins lonal  array  FM;  note  that  the 
order  of  integration  and  summation  has  been  reversed  for 
ease  in  computation.  The  lopic  after  statement  number  40 
is  used  to  compute  the  phase  error  variance  for  QCM  if  a 

2 

fla^r  (IQC^>0)  is  declared.  The  averaf^e  value  of  I and  I 

is  then  obtained  for  all  possibilities  within  the  coherence 

I 

' time  of  the  loop  (for  this  thesis  equivalent  to  ten  data 

pulses',  and  the  average  phase  error  variance  is  computed 
I for  500  s i/rnal -enerpy-to-no ise  ratio  points. 


j 

I 

I 

I 

i 

t 

I 


I 
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Table  XII,  continued 
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(20.  cont.)  ■ 

carrier  tracking  capability  of  a QCM  scheme  is  equal  to  that  of 
BPSK  for  the  same  closed-loop  bandwidth  and  time  bandwidth  prod- 
uct (BT)  of  the  channel  filter  at  a specified  signal  energy  to 
noise  ratio  (Sb/No). 

The  data  detection  performance  for  QCM,  BPSK,  and  QPSK  is 
analyzed  by  comparing  the  one-shot  probability  of  error  con- 
ditioned on  a phase  error  as  a function  of  Eb/i'Jo  and  BT  for 
a specified  channel  filter.  Carrier  tracking  performance  for 
QCM  and  BPSK  is  analyzed  by  obtaining  an  average  phase  error 
variance  for  the  linear  model  of  a Costas  loop.  In  computing 
both  the  probability  of  error  and  phase  error  variance,  the 
intersymbol  interference  is  modeled  from  a truncated  data  se- 
quence. In  addition,  the  bit  energy  penalty  for  tracking  a 
QPSK  signal  versus  a BPSK  or  QCM  signal  with  no  ISI  is  examined. 
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